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1 Introduction 

In recent years there has been considerable progress in extending the holographic corre- 
spondence to spacetimes that are not asymptotically AdS. Schroedinger [1, 2] and Lifshitz 
[3, 4] spacetimes were introduced as possible duals of non-relativistic strongly coupled field 
theories. 

Holographic renormalization for these spacetimes was explored in [5] for the Schroding- 
er case and more recently in [6-8] for the Lifshitz case. In particular, it was suggested in 
[7] that the non-relativistic analog of the Weyl anomaly, derived in the context of ordinary 
AdS/CFT in [9], might be present for 3-1-1 dimensional Lifshitz spacetimes with z = 2. 
Conformal anomalies play an important role in relativistic field theories, especially in two 
dimensions, where various physical quantities display universal behavior which is governed 
by the central charge only. It is clearly of interest to explore to what extent similar results 
carry over to non-relativistic field theories. 
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Motivated by these observations, in this paper we explore the general structure of the 
Weyl anomaly for non-relativistic field theories in d = 2 + 1 dimensions with z = 2. The 
cases d = 4, z = 3 and d = 6 were analyzed in [10] and [11] respectively. Contrary to the 
relativistic case, where anomalies are present only for even dimension, we show that in the 
non-relativistic setting anomalies can be generated in odd dimension as well. That this is 
in principle possible can be quite easily seen by dimensional analysis. The analogue of the 
trace of the stress-tensor for non-relativistic field theories has dimension z + d— 1. A term 
with a time derivatives and h spatial derivatives, on the other hand, has dimension az + h. 
For generic z, there can only be contributions to the conformal anomaly with a = 1 and 
b = d — 1. Such terms can only appear in theories which explicitly break time reversal 
invariance, which will not be the case for the theories that we consider. For special values 
of z, especially for integer values of z, other values of a, 6 are allowed: for example, for 
d = z + 1, terms with either (a, 6) = (2,0) or {a,b) = (0,22:) can appear. In particular, 
for the case of interest to us with z = 2, d = 3, one can have terms with either two time 
or four spatial derivatives. There are many such terms that one can write down, but the 
more detailed analysis that we give below shows that the anomaly is generated by a total 
of two linearly independent and non-trivial structures. 

In order to show that the anomaly is indeed generated in field theory models enjoying 
Lifshitz symmetry, we compute the heat-kernel expansion of a fourth-order differential 
operator that transforms covariantly under local anisotropic rescalings. The corresponding 
free field theory model has various applications in condensed matter physics, see [12] and 
references therein. We explicitly compute the numerical coefficients that appear in front of 
the two possible contributions to the anomaly, and find that only one of them is nonzero. 

An obvious question is then whether a similar statement holds for those purported 
strongly coupled field theories which are dual to Lifshitz spacetimes. We compute the two 
central charges using the holographic renormalization methods of [7]. Just like in the free 
field theory example, only one of the coefficients turns out to be non- vanishing. The results 
for the holographic computation perfectly agree with those of the recent paper [13], which 
computes the full conformal anomaly using the methods of [6] and which appeared while 
we were preparing this paper for submission to the ArXiv. 

Remarkably, the anomaly in both models is entirely generated by terms that only 
contain time derivatives, even though more general structures containing space derivatives 
are in principle allowed. It is not clear to us whether there is a simple explanation of this 
fact. 

The outline of this paper is as follows. In section 2, we setup and perform the field 
theory computation of the conformal anomaly. The corresponding analysis for Lifshitz 
spacetimes is presented in section 3, and in section 4 we present our conclusion. Some 
technical details are discussed in the appendices, such as the complete classification of the 
possible structures that can appear in the conformal anomaly, and an interpretation of the 
bulk metric and gauge field in the field theory language. 

Note added: In an earlier version of this paper, we erroneously claimed that there 
were three instead of two possible conformal anomalies and used holographic renormaliza- 
tion to compute two of the three central charges. The aforementioned paper [13] correctly 
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pointed out that there can only be two and not three contributions to the anomaly. We 
corrected this mistake in this version and as a result our original anomaly computation 
now yields the complete anomaly. 

2 Anisotropic-Scaling Anomaly in the Quantum Lifshitz Model 

Scale invariance arises naturally when describing systems at a quantum critical point. Even 
when the underlying theory is not Lorentz invariant, it may still exhibit invariance under 
so-called Lifshitz-like scaling that treats time and space differently, 

x^Xx, t^yt. (2.1) 

Only the specific value z = 1 is compatible with Lorentz symmetry. In the high-energy- 
physics community (2.1) is often called anisotropic scaling, though in the consensed-matter 
literature the word 'anisotropy' is typically used for spatial anisotropy only. 

In this section we analyze the anomalous breaking of anisotropic scaling symmetry in 
a three-dimensional field theory known as the quantum Lifshitz model [12]. This model 
arises in the classical field theory description of three-dimensional Lifshitz points, which is 
why the corresponding quantized theory is called the 'quantum Lifshitz model'. We refer 
the interested reader to [12] and references therein for a discussion on applications of the 
quantum Lifshitz model to condensed matter systems such as the so-called quantum dimer 
model. 

We shall find that the anomalous breaking of the Lifshitz-like scaling (2.1), which we 
will sometimes refer to as the 'Lifshitz anomaly', is parametrized by two central charges. 
In this section we will use heat kernel methods in order to determine these two central 
charges in the quantum Lifshitz model. 

2.1 The Quantum Lifshitz model 

We consider a three dimensional free boson (j){t, x) with the following Euclidean action 

S = j dtd^xNVh (^^N'\dt(t>)^ + \i^(t>?^ ' (2-2) 

where we have coupled the system to the degrees of freedom of an auxiliary three-dimensional 
metric with a preferred time foliation, of the form:-^ 

ds2 ^ jv2^^2 _^ hijdx'dx^. (2.3) 

Furthermore, A = VjV* is constructed out of the covariant derivatives of the spatial metric 
hij. The metric (2.3) keeps its form under diffeomorphisms in time, 1 1— )■ T(i), and in space, 
I— 7- i^*(x), and both leave the action (2.2) invariant if we transform N and hij as dictated 
by (2.1). This is one of the reasons why it is convenient to combine N and hij in an auxiliary 

^We did not include a shift A^' , because it can be removed locally by a foliation preserving diffeomor- 
phism, e.g. 1 1— >■ r(t) and t-^ ^'(t, £). 
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metric, the other reason is that this structure also naturally emerges in the holographic 
setup. 

By integrating by parts and ignoring boundary terms^, the action can be written as: 

S = j dt(fxN^/h(l)D(p, (2.4) 

where D is given by: 

D = -^^dtN~^Vhdt + ^ANA (2.5) 
This model is classically invariant under local anisotropic scale transformations, that is 

^ e^'^N, hij e^'^hij, (p ^ (p, (2.6) 

where w is an arbitrary function of t and rr*. In particular, the operator D transforms as 

D e-^'^D (2.7) 
The variation of the action S under a Weyl transformation is thus given by 

SS = / .M^. [m^ . = / ,t,^.S, [2N§ + 2.„i| 

Defining the energy density £ and momentum flux (spatial stress tensor) Iljj as 

we see that local anisotropic scale invariance implies: 

2f + = 0, (2.10) 

which is the non-relativistic analog of the tracelessness condition = 0. 

We will prove in this section that the classical anisotropic conformal invariance of this 
model is broken at the quantum level. The quantum expectation values of the energy 
density {£) and the spatial stress tensor (|n*-'^ are given by (2.9) when one replaces S by 
the effective action W = — log Z. In the presence of an anomaly, the right-hand side of 

2{£) + (jV^=A (2.11) 

is non-zero, where A is the anomaly. Equivalently, one may express the anomaly directly 
in terms of the variation of the (renormalized) effective action W , 

SW = j N^/h6pA, (2.12) 

where 5p = uj is the infinitesimal scale factor. We will show that 

A = (iGh'mdtiN-'hij) + 5{h'^hi,f - lOh'^h.kh'^'hii^ + ^lviJ^ (2.13) 

^In the following, we will assume that the theory is defined on a manifold without boundary. 
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where a dot indicates dt and VjJ* is a "trivial" total derivative, by which we mean that 
it can be removed by adding appropriate local counterterms. In fact we will show that by 
adding appropriate counterterms to the action (2.2), the anomaly can be written in the 
simpler form 

This anomaly will be computed using a heat-kernel expansion. 

As an aside, one of the reasons for the particular interest in this model is that the 
ground-state wave functional is invariant under time-independent conformal transforma- 
tions in space. All equal-time correlators can be computed using the machinery of a 
two-dimensional field theory [12, 14]. One may thus naively expect that the anomalous 
breaking of anisotropic scaling symmetry (2.1) is somehow related to the two-dimensional 
Weyl anomaly (T*j) oc R, see e.g. §5. A of [15]. We find, however, that this is not the case. 
The anomaly that we find involves only derivatives with respect to the time coordinate, 
whereas the two-dimensional Ricci scalar R obviously only contains spatial derivatives. 

2.2 Heat-kernel expansion 

The quantum effective action W for the model (2.2) can be computed explicitly, and is 
given by the formal expression: 

W (x^lndet{D), (2.15) 

where det(I?) is the determinant of the operator D defined in equation (2.5). As usual, 
this determinant is not well-defined and must be regularized. We will employ (^-function 
regularization. We define the generalized zeta function as 

C{sJ,D) = TrL2ifD'^), (2.16) 

where s is an arbitrary positive number and an appropriate function space on which 
is trace-class. The regularized effective action is given by [16]: 

W = -ic'(0, hD)-^ Hf^'nO, 1, D), (2.17) 

where C'(0) /> = dsCis, f, D)\s=o and /i is the usual arbitrary renormalization scale. The 
zeta function (^(s,f,D) is related via a Mellin transformation to the heat kernel:"^ 

K(e,/,Z5) = Tri2(/e-^^), (2.18) 

where / is an arbitrary function of t and x^, and e is an arbitrary positive parameter. In 
principle K depends on the global behavior of the operator D (the trace can be written 

^In particular, the relation between the two is 

as,f,D) = r{s)-' 1^°^ dee'-'K{e,f,D), K{e,f,D) = ±. <ji ds e'^ r{s) ({s, f, D). 
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as a sum over the spectrum of the operator, which is determined by global properties); 
however there is an asymptotic series of the form: 

oo 

Kie,f,D)^^e'2-^dk{f,D), (2.19) 

A;=0 

where dk{f,D) can be computed locally from N and hij. By repeating the analysis of 
[16] section 7.1, it is possible to show that the variation of the renormalized effective 
action under an infinitesimal anisotropic local scale transformation /i — )■ (1 + 25p)h, N — )• 
(1 + 2Sp)N, is given by ^: 

6W = -2d2{6p,D). (2.20) 

In other words, the anomaly is given by the term in the heat-kernel expansion. As 
explained above, this will be a local functional of N and h; we will therefore write: 

d2{f,D)= I dt(fxNVhfa2{N,hij), (2.21) 



where 02 (A^, hij) is a local function that depends on N and hij. The Weyl variation of the 
marginal heat kernel coefficient 02(1, -D) vanishes. To see why, consider the Weyl variation 
of the heat kernel coefficient at order k, 



d 

dj 



dkihe-^-'^D) = i4-2k)dk{f,D), (2.22) 

7=0 



which vanishes when k = 2. We used the identity for the Weyl variation of the full heat 
kernel, 



d 

d'y 



K{e,l,e-'^fD)= j- 
7=0 ^7 



r|.^/g-ecxp(-47/)D\ (2.23) 

7=0 

= 4eTr(/L>e"'-°) (2.24) 

= -46^Tr(/e-^^) (2.25) 

= -AejK{eJ,D). (2.26) 

As a consequence, we can use this property that 5a2(l,-D) = to find that, in local 
language, 

NVha2{N,hij) = e^"^ A^\//Ia2(e^'^ A, e^'^/ii^) + total derivatives. (2.27) 

In other words, the anomaly itself must be Weyl invariant (up to total derivatives), and 
this puts very strong restrictions on the possible terms that can appear in 02. 



*The factor 2 comes from the factor 4 in D e '^''D under scale transformations. 
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2.3 Analysis of the possible terms 

Just by dimensional analysis (i.e. by requiring invariance under constant rescalings) we can 
see that there are many terms of the right dimension that can appear in the heat kernel 
expansion. In dimensional analysis, dt has dimension two and di has dimension one; we 
are interested in terms of dimension four that are covariant under reparametrizations of t 
and reparametrizations of x*. These reparametrizations should not mix t and since that 
would ruin the form of D. Examples of terms of the right dimension are^ 

N-^dthijG'^'^'dthki, 

R^, N-'^ANR, (Af-iAiV)2, 

(2.28) 

h'^h''^{N~^diN){N~^dkN){N~^VjdiN). 



where G^^^'- = \{g^^g^^ + g^^g^^) — Xg'^^ g^^ is the DeWitt metric and A an arbitrary real 
number (in General Relativity A = jy^-, where D is the spacetime dimension). Terms 
with one time derivative and two space derivatives cannot appear because of time reversal 
symmetry. All these terms are scale invariant, but the anomaly should also obey the Wess- 
Zumino consistency condition (namely the second variation of the effective action W should 
be symmetric). The possible consistent structures are classified in Appendix A. The most 
general form of the anomaly is therefore: 



Ian + 0(^8.n) (1 



+ C2 [R- —AN + h'^ ( —diN ] ( —djN ] ] (2.29) 



+ (trivial total derivatives). 
Here the trivial total derivatives can be cancelled by appropriate counterterms. 
2.4 Computation of the anomaly 

As we explained above, the spectral function can be expanded as: 

K{e, f,D) = Y^ et"i f dt(fxNVhf{t, x)ak{N, hi^), (2.30) 
fc>o 

where ak{N,hij) is a local function of N and hij. To evaluate this we need a suitable 
basis; it is customary to use the rescaled Fourier modes so that they are orthonormal with 
respect to the measure that includes the Ny/h factor. Nevertheless, as pointed out in [17], 
the cyclicity of the trace allows us to use the usual flat Fourier modes. We thus find 

K = J J did2^e-^"*-''="/e-^^e*'^*+*'=^ (2.31) 



complete classification can be found in Appendix A. 



-7- 



We can conjugate the Fourier mode to the left to get the expression 

du}(Pk 



K 



j^j "td'^fe-"". (2.32) 



where D2 is obtained from D by shifting the derivatives as follows: 

dt ^ dt + ioj di ^ di + iki. (2.33) 

The most singular term in the heat kernel is the one where we keep only the terms in D2 
without derivatives, leading to 

'MLD) = f^f dtd^xfe-<^-'-'^^''^'\ (2.34) 



e J [zn 

where k'^ = h^^kikj. This expression is readily evaluated to yield the first term in the heat 
kernel expansion: 

do{f,D) = I dtd^xNVhf{t,x) (2.35) 
167r J 

Computing the subleading terms in the heat kernel expansion is now straightforward though 
somewhat involved. We shall write 

D2 = Dl + Dlf-^ (2.36) 

where D2 is the piece we isolated above that contains and k^, and D!f^ the remainder. 
We then expand the exponential of -02°*- contains a factor of e, but u counts as e~^/^ 
and k as e~^^^ in the Gaussian integral, therefore -D™* has a term which scales as e~^/^, and 
to get to the finite term one needs to expand D™* up to fourth order so that we get terms 
up to k^"^. However, the problem becomes tractable if we consider the time-derivative and 
space-derivative sectors separately. This is consistent because the anomaly can only have 
structures involving either two time derivatives or four spatial derivatives. 

The two-derivative anomaly 

In order to compute the two-derivative contribution to the anomaly, and in turn Ci, it is 
sufficient to consider metrics that only depend on t, and not on x"". Thus we can drop all 
the terms with spatial derivatives di in Z)™*. Moreover, by changing the coordinate t if 
necessary, we can take = 1. With these assumptions, we have: 

D^ = u^ + {k^f, (2.37) 

Df = -iuj{dt + -^dtVh) - ^dtVhdt. (2.38) 
yh yh 

We have to expand to second order in D™*- Since 1)2^* and D2 do not commute, we use 
the following formula: 



JQ<a<l Jo<a+B<l 



I0<a+I3<l 

(2.39) 



-8- 



We find the following contribution to a2-^ 

a2{f, D) = j dtcfxVhf{l6h'^hij + m'kjf - lQh'^hjuh'''hi,] + ... (2.40) 

where the ellipses denote possible four-derivative contributions. To reinstate A^, we simply 
need to change dt — )• dtN and dt — )• N~^dt- We can remove the first term by adding local 
counterterms as explained in Appendix A, to obtain the two-derivative contribution to the 
anomaly:^ 

^—4? (h'^h''^ hikhii - -{h'^hi.)A . (2.41) 
1287riV2 \^ j< 2^ y ^ ^ 

Using (2.20) and (2.29), we see that 

Ci = ^— . (2.42) 
1287r ^ ^ 

The four-derivative anomaly 

We now determine the four derivative contribution, and in turn C2. As explained in Ap- 
pendix A, there are 6 possible terms that can appear, 5 of which are total derivatives. These 
structures are distinguished by a metric of the form hij = e^^^^^Sij and N = e^^^\ which 
can be used to simplify considerably the computation. The four-derivative contribution to 
the anomaly is thus 



^ = 4^^^^ (-5(5^iV)i? + 3(5^iV)(lAiV) + 2{d^N){^V,d^N) - 5d^AN ) . (2.43) 



It is interesting to note that this result is a total derivative, and as predicted by the Wess- 
Zumino consistency condition, it is orthogonal* to the non-trivial total derivative defined 
in equation (A. 15). As a consequence, this term can be removed by a local counterterm 
and we conclude that 

C2 = 0. (2.44) 
In Appendix C we present an alternative derivation of C2 = 0. 

The anomaly 

In summary, the Lifshitz model (2.2) exhibits an anomaly under anisotropic local scale 
transformations, which after the addition of appropriate counterterms is given by: 

^ = {h^'h'^' K,h,i - \{h^^h,f^ . (2.45) 

It is striking that the anomaly involves only time derivatives. So far, it is unclear to us why 
this happens. It is also in contrast to the naive expectation that the anomaly is somehow 
related to the trace anomaly of a two-dimensional conformal field theory, as we mentioned 
at the beginning of this section. 



®This computation is in principle quite lengthy. However, since there are only few terms that can appear, 
one can work this out for a diagonal hij and then reconstruct the full answer. 
''There is a relative factor -2, see (2.20). 

*To see what we mean by 'orthogonal', we refer to the appendix A. 
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3 Holographic Calculation 



In the previous section we showed that a theory with anisotropic scahng symmetry has an 
anisotropic scahng symmetry anomaly parametrized by two central charges, denoted by Ci 
and C2. We computed these central charges for a particular model defined by the action 
(2.2). In this section we show that these central charges can be computed holographically 
for the Lifshitz spacetime considered in [6, 7]. 

3.1 The Hamilton Jacobi Formalism 

Lifshitz spacetime is a proposed gravitational dual to a field theory at a UV fixed point 
with anisotropic (Lifshitz-like) scaling symmetry t 1— )• A^t, x 1— )■ Ax. The configuration of 
{d + l)-dimensional Lifshitz spacetime that we consider is a solution of the field equations 
derived from the Einstein-Proca action, 



(3.1) 



The parameters of the theory are the curvature scale i and the dynamical exponent z; they 
are related to the cosmological constant and the Proca mass via 

_ z2 + (d-2)z + (ci-l)^ 2 _ id-l)z 

We shall work in units such that IGvrG = 1 and i = 1 henceforth. Also, the computation 
that we shall perform in this section shall be for the case of d = 3 boundary dimensions, 
though we would like to keep d explicit where this may be illuminating. The solution is 
then given by 

ds^ = dr"^ -e^^'df + e^''dx^, A = V^e'''dt. (3.3) 

with ao = —{z — 1) jz. A shift in the radial coordinate r 1— >• r + log A generates the Lifshitz 
scaling t ^ X^t, x ^ Ax. 

The goal of this section is to compute the (divergent piece of the) on-shell value of the 
above action. We shall do this in the Hamilton-Jacobi (HJ) formalism [18, 19]. Much of this 
section should be seen as an extension to [7], so we refer to that work for a more detailed 
discussion of the HJ formalism. For instance, the ADM-like Hamiltonian associated to 
(3.1) was derived in [7]. The HJ equation is a differential equation for the on-shell action, 
i.e. solving the HJ equation will give us the on-shell action 5ci. We write the HJ equation 
H = Oas^ 

{5ci,5ci}-£ = (3.4) 



^Again, see [7] for more details. In particular, (3.4) is the Hamiltonian constraint; the momentum 
constraint is automatically satisfied by choosing an Ansatz for (the local part of) Sd that is covariant on 
the radial cut-off slice. 
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where the brackets are given by the fohowing expression. Let F and G be two arbitrary 
phase-space functionals, then the brackets are defined as 



{F,G} = 



1 \ 6F SG I 6F 6G . . 

laclbd - -1 Tlablcd T T h l^lahYlT-FJ' ^^-^l 

d-1 ) b-fab hcd 2 bAa OAi, 

J_ 6F ^ 6G 



where 'jab and Aa are the induced fields puhed back onto the radial cut-off slice. The 
brackets (3.5) were introduced in [18]; see also [7]. It should be noted that these brackets 
are only introduced as a short-hand notation for the 'kinetic' part of the Hamiltonian 
contraint; they are not Poisson brackets (or any other type of special brackets). 

We split up the on-shell action Sd into a local power-law divergent piece 5'ioc and a 
non-local piece W that is at most logarithmically divergent; W is finite in the absence of 
anomalies. Using the split = jSioc + W, we may write the HJ equation as 

= {Sloe, Sloe} -C + 2{Sioc,W} + {W,W}. (3.6) 

We define Tiioc = {Sioc, Sioc} — ^, and we require that the divergent terms cancel.^'^ It might 
happen that we get a finite remainder Urcm that cannot be removed by local counterterms. 
As a consequence, the HJ equations in the large r limit, where we can ignore {W, W}, give: 

2{5ioc, W} « -Tirem. (3.7) 

The symbol "~" means an equality in the large r limit. We will presently see that this is 
related to the anisotropic scaling anomaly in the dual non-relativistic field theory. First, 
we need to specify the boundary conditions in order to be able to identify which terms are 
divergent. 

3.2 Boundary conditions and anomaly 

From the field theory side, we know that the volume form has a definite scaling weight 
[Vol(i] = [dtd'^~^x] = z + d — I. In the dual gravitational picture this weight is translated 
to a radial scaling, such that 

z + d-1 = [Volrf] = = drhgN + -drlogdeth (3.8) 

Here, N is the lapse function and hij is the induced metric on a spatial slice of S^. We 
assume that the spatial metric is of the form hij = e^^hij, where hij has a finite limit as 
r — )■ oo. It then follows that the lapse scales as ~ e^*". This puts a restriction on the 
degrees of freedom contained in the metric. It implies in particular that we must turn off 
the off-diagonal mode in ^ti that scales as e^^''; in terms of the linearized modes discussed 
in [6, 7, 20], one needs to kill the cu mode. This naturally leads us to consider only 
deformations with a preferred time foliation, as suggested also in [6]. 



Here we refer to a term F as "divergent" if limr->cx) ^/^F — oo. 
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Let us redefine N and hij to be the renormalized lapse and induced metric/^ 

N = lim e-^''(-7**)-^/2^ (3.9) 

r—>-oc 

hij = lim e"^''7jj, (3.10) 



¥00 



such that the renormalized volume form is given by NVh = lim^^oo e-'^'+'^-'^> With 
these conditions, it is straightforward to see that: 

n6W 6W - 6W 

2^1 {Sloe, W} « 2ziV2^ + ^hij— + zAtjj « - V7^rcm, (3.11) 

where At = limr^oo^~^^At- As noted in [6, 7], for z > 2 and d = 3 the vector field 
becomes an irrelevant operator, which prevents us from defining the on-shell action non- 
perturbatively in the sources. In particular, the linearized analysis for z = 2 shows the 
presence of a mode that shifts the value of a = A'^Aa and introduces logarithmic divergences 
in the metric: 

Itt = -e^' {l + (4c3 ^ + C4) + ^ (4ci - 5C2 + 4c2 r) e"^'^^ , (3.12) 
m = (i^'^^ii ( 1 + (-2c3 ^ + C5) + ^ (4ci + 5c2 + 4c2 r) e"^^ 



a = -1 + 2c3 + 2 (ci + C2 r) e 



-4r 



It is convenient to expand the counterterms in powers of a + 1, i.e. in fluctuations of 
a around its background value ao = —1. However, turning on the mode C3 means that 
(a— ao)"" is not suppressed as r — )• 00, and an infinite set of counterterms would be required. 
Furthermore, the mode C3 introduces logarithmic divergences in the metric sector, spoiling 
our definition of anisotropic conformal infinity (3.9). For this reason, we will turn off this 
mode^^ by setting C3 = 0. In particular notice the important relation: 

At = e^^'N + subleading, (3.13) 

where the subleading terms scale as e"^''. In the frame field language of [6], this corresponds 
to 5{Aa) = 6{Aae'\) = 0. Notice that with these boundary conditions, W becomes a 
functional of N and h only. Therefore we have: 



dW 



dN 



_ 6W_ 6W dAt 



where the variations on the right are unconstrained, while the variation on the left repre- 
sents the total variation of W with respect to A^. Therefore (3.11) becomes: 

dW 5W , 

2^7^ + = -\^^rem. (3.15) 



^^Again, we do not consider a shift A*'*, as it can locally be removed by a foliation-preserving diffeomor- 
phism. 

^^As shown in [7], it is possible to consistently impose this condition when higher order non- linear 
corrections are considered, and we believe there is no obstruction at the full non-linear level. 
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By comparing with (2.11), we find that the anomaly is given by 

A= - Wui e^^'Urcra- (3.16) 

r— >oo 

This is the holographic anomaly in the HJ formalism. In the following we will use the 
aforementioned boundary conditions in order to identify possible divergent terms. 

3.3 Two-derivative counterterms and anomaly 

The Hamiltonian constraint at the level of two spacetime derivatives is given by 

= {"^loc, Sioc}^') -R+ (3.17) 

Now, let us take the most general Ansatz for the on-shell action in such a way that T-La = 
(i.e. covariant on S,.). At the non-derivative level, we thus take 

s[oc = I d'^x^Foia), (3.18) 

with FQ{a) an arbitrary function of the Lorentz scalar a = AaA'^. This case was analyzed 
in detail in [7], and does not lead to anomalies in the present case. At the level of two 
spacetime derivatives, we take^^ 

^loc = / d^'x^ {Fi R + F2 DaAh D'^A'' + F3 DaA^ D^'A'^ + J4 {DaA'^f + F^ d'^a daa 
+ Fq A^DhAa A'^D.A" + F7 A" daa Df,A^ + Fs A^DaA^ d^a + F^ {A'^daaf^, 

(3.19) 

where each coefficient Fi is a function of a. We do not include terms that involve sec- 
ond derivatives, e.g. ~ DaDj,Ac, since those can be expressed in terms of first-derivative 
terms by means of partial integration. Equation (3.17) can be solved straightforwardly by 
expanding the coefficient-functions around the Lifshitz background, 

-^^(«) = Yl /^(n)(«-«o)". (3.20) 

n>0 

(2) 

The coefficients are then found by solving the equation TiU^ = recursively order by 
order. 

For the case of 2; = 2 in d = 3 boundary dimensions one finds a break-down of the 

(2) 

recursive 'descent' equations, i.e. ^-1^^ 7^ for any choice of coefficient functions F^. Such 
a break-down signals the presence of an anomalous remainder Tirem giving rise to the 
holographic anomaly. As explained in the introduction to this section, we should still be 

^^In order to get a general understanding of these terms, notice that the terms proportional to Fi with 
i = 1, 2, 3 are of the schematic form DA DA, for i = 4, 5, 6, 7 they look like AA-y DA DA, and for i — 7 
is AAAADADA. 
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able to determine the coefficient of the divergent part of the on-shell action in order to 
have a weh-defined theory on the boundary. Doing so leads to the following remainder i^'^ 

nill = -DaA.D'A^ + ^iDaA'^f, (3.21) 

which is obviously affected by the ambiguity of adding finite local counterterms to the 
action. Using the boundary conditions described in 3.2 and the definition of the holographic 
anomaly (3.16), we find the following two-derivative contribution to the anomaly: 

^ = ^ ^2 (^*'^'' h.'^f'^l - + • • • > (3-22) 

where the ellipses denote possible four-derivative contributions. The above anomaly is 
indeed of the form (2.29) with: 

We reinstated the four-dimensional Newton's constant G and the curvature length scale i. 
Also, hij is generically related to the extrinsic curvature by Kij = {hij — ViNj — VjNi)/2N 
(and we had set the shift iVj to zero). 



3.4 Four-derivative anomaly 

One may repeat the above steps at the level of four derivatives.^^ The four-derivative 
Ansatz is 

^loc = / ^''^Vi {Si + Q2 RabR"' + G3 RabcdR"'"' + GaOr) + ■ ■ ■ (3.24) 

The ellipses denote terms that involve the Proca field Aa- All the terms that appear at 
this level are finite for our choice of boundary conditions, which means that they can only 
contribute with trivial total derivatives to the anomaly. In this case, we find the remainder 

^^e^ = lR'-\Ral>R''' + ... (3.25) 

where the ellipses denote once again terms that involve Aa- Writing this in terms of the 
two-dimensional Ricci tensor gives: 

nitl = l%,'R''-leRf + ..., (3.26) 

where these and '^Rij are the too-dimensional Ricci scalar and tensor and we have not 
written down terms that involve derivatives acting on N. We can use the off-shell identity 

To see how this comes about, one can can look at the solution Sioc = / d'^a;-^— 7 /Iioc for values of the 
dynamical exponent close to the critical value, 2 ~ 2. One finds a simple pole at 2 = 2, 

A*o^ = j^{-DaA,D'A- + ^{D^A''f)+0{{z-2f), 

where the residue is simply the two-derivative remainder (3.21). 

^^One does not expect to find anomalous contributions that contain three derivatives (one time and two 
spatial), since terms that involve an odd number of time-derivatives are not invariant under time-reversal. 
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that relates the Ricci tensor to the Ricci scalar, = ^"^Rhij, which is specific to two 
dimensions. When we plug this into (3.26), we find that the Ricci-squared terms cancel 
and we have (equation (2.29)): 

C2 = 0, (3.27) 

which agrees with the field theory computation. 

Notice that while we were able to extract the coefficient C2, we have not performed a 
complete analysis of the counterterms at the four-derivative level, which would be rather 
involved. Nevertheless, the complete answer has been computed using the results of [6] in 
[13], in perfect agreement with our result C2 = 0. In conclusion, the holographic anomaly 
is given by 

-4 = ^ ^ {h^'h'^' h^kkl - \{^'^h,,f) ■ (3.28) 
4 Discussion and conclusions 

In this paper we computed the anisotropic scaling anomaly of two Lifshitz theories, one 
defined using a standard field theory quantization of an explicit classical action (2.2), 
the other defined using the holographic correspondence. A precise definition of Lifshitz 
holography is still lacking, and a microscopic definition of the strongly coupled field theory 
is not known. It is therefore a priori not very meaningful to compare the two anomalies. 
Nevertheless, we found that the anomalies are quite similar. In both cases there are two 
possible central charges of which one vanishes, and as a consequence the two anomalies are 
directly proportional to each other. The ratio of the two anomalies, in the conventions used 
in this paper, is 21/ G, with I the curvature radius of the Lifshitz spacetime and G the 4d 
Newton constant. It would be interesting to evaluate this quantity in explicit string theory 
embeddings of Lifshitz spacetimes to see how it scales with the various integer fluxes, as 
this will provide some measure of the effective number of degrees of freedom of the dual 
field theory. 

It is quite mysterious that the conformal anomaly only involves time derivatives, it 
is even mysterious that there exists a conformal anomaly at all. According to [21], the 
dynamical critical exponent is in general renormalized, and as soon as 2; = 2 + e a conformal 
anomaly can no longer be written down. So either there is some unknown mechanism that 
protects the value of z = 2, or the conformal anomaly can be removed in the full quantum 
theory. In the latter case, one would be in the peculiar situation that one would need to 
include counterterms that diverge in the classical limit. Further work will be required to 
clarify this issue. 

It is also clearly of interest to explore other systems with anisotropic scale invariance 
to examine whether the conformal anomaly is still of the same form. In particular, when- 
ever one has a Lifshitz solution in a theory with Chern-Simons type terms, time reversal 
symmetry is broken and it is logically possible to have contributions with an odd number 
of time derivatives to the conformal anomaly. It is in principle straightforward to extend 
the analysis in appendix A to determine whether there are non-trivial terms of this type 
and we leave this as an exercise. 
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As mentioned before, one of the main uses of the conformal anomaly is that it is a 
relatively simple property of a field theory which sometimes gives rise to certain universal 
properties. For example, in the relativistic case, in d = 2 the conformal anomaly com- 
pletely fixes the free energy at high temperatures, and the central charges also control the 
logarithmic divergences in the entanglement entropy in d = 2, 4. Whether similar universal 
properties also exist for non-relativistic field theories is an interesting open problem that 
we hope to come back to in the future. 
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A Classification of possible terms in the anomaly 

In this appendix we explore to what extent it is possible to remove total derivatives from 
the anomaly. This is achieved by adding appropriate scale invariant counterterms to the 
action that are not invariant under local scale transformations. Clearly, we can discuss 
the two-derivative and the four-derivative terms separately. Let us start with the former; 
there are only three possible scale-invariant terms that we can construct with two time 
derivatives: 

h^'j^dti-dth,,), —(h'^h,,f, h^'h.kh'^'hi,. (A.l) 
It is straightforward to see that the two combinations 

h'^]^dt{^dthii), ^ [h'^h'^^ Kuh^i - \{h'^hjf^ , (A.2) 

are invariant under local scale transformations (up to total derivatives). These two terms 
are related by partial integration, and we now show that it is indeed possible to "partially 
integrate" inside the anomaly by adding an appropriate counterterm to the action. The 
most general form of the anomaly at the two derivative level is: 

5W = j 5p^a^h'^dt{^hij) + [h'^h^' hikh.i - ^{h'^k.f^ I . (A.3) 

The presence of the factor 5p prevents us from doing partial integration directly. Let us 
add the following counterterm to the action: 

W' = W + c [ NVh^{h}^kjf. (A.4) 



-16- 



It is then easy to check that 

5W' = j 5p\{a- 8c)^h'^dt^k, + {b + 8c)^ (h'^h^' hikh.i - ^{h'^k.f^ | . (A.5) 

Therefore we can pick c = a/8 and get rid of the first term, which is tantamount to 
integrating by parts, or discarding total derivatives in the anomaly. For instance, in the 
field theory analysis, we went from (2.40) to (2.41) using this procedure. In particular, we 
had a = l/487r and b = -5/3847r, such that b + 8c = a + b = I/UStt}^ 

Let us now consider the four derivative level. In this case we are interested in terms 
of the form Vj J* in the anomaly. We ask ourselves to what extent it is possible to remove 
them by adding local counterterms G to the action. Both the total derivatives and the 
local counterterms must be scale invariant, therefore there is only a finite number of them. 
Let us choose a basis: 

4 a = l,...,N (A.6) 
Gb b = l,...,M. (A.7) 

The Weyl variation of a linear combination qh Gb can be written, after partial integra- 
tion, as: 

<5j;%G6 = ^^Af,bgbV,Jl, (A.8) 

b ab 

If the variation of the effective action reads: 

SW = J Vhuj + CaVi Jl j , (A.9) 

we can get rid of the total derivatives if we can solve the system of linear equations: 

Mabqb = Ca. (A.IO) 

If we are to remove all the possible total derivatives that can appear, the number of rows 
N of the matrix Mab must be less than or equal to the number of columns M, and the rank 
of the matrix should be maximal. It is easy to check that there are 6 possible functionally 
independent scale invariant currents J*, and we choose the following basis: 



Ji = Nd'R 4 = {d'N)R 4 = {d'N){j^djN){j^dm) 

4 = {d^N){^AN) 4 = {d^N){jjVjd^N) 4 = d^AN. 



^^Bear in mind that there was an extra factor —2 coming from the relation between 0.2(3 p,D) and the 
integrated anomaly, cf. (2.20). 
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Analogously, there are 12 functionally independent scale invariant counterterms, and we 
choose the basis: 



G4 = (ia,iv)(iaw)i? 



G2 = AR G3 = {jrAN)R 

G, = (ij^d,N){^d^N)f G, = ij^d,N){^d'N){j^AN) 



G 



10 



±v,d,NijV'dm Gu = -hA'N 



N 



G 



12 



(A.12) 

While we have many more possible counterterms than currents, it is important to stress 
that not all the counterterms are independent, since we can always partially integrate 
inside the action. This means that some linear combinations of counterterms will have 
the same Weyl transformation. Furthermore, there can be Weyl invariant combinations of 
counterterms that do not help in removing total derivatives from the anomaly. 

By taking the Weyl variation of the 12 terms Gb, it is straightforward to compute the 
matrix M, which is given by: 



/-4 


2 


2 
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-4 


-2 


-2 


-4 
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-8 
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-5 


-6 
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-8 
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-2 














-4 
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-2 
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-4 
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-2 


-2 











4 





-4 
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2 / 



(A.13) 



It is easily checked that M does not have maximal rank (which would be 6), but it has 
rank 5. In fact, M has a 7 dimensional space of null vectors, which is spanned by the 6 
total derivatives Vj J* and a Weyl invariant term: 



6 I VhViJ'a 



0, 



NVh{R+^AN-j, 



diNd'N 



0. 



(A.14) 



Since the rank of M is 5, the Weyl variation of the most general counterterm spans a 5 
dimensional subspace of the 6 dimensional space generated by VjJ*. That means that 
we can find an orthonormal basis (with respect to the usual Euclidean scalar product 6ab) 
for the currents where 5 are trivial (i.e. removable by counterterms) and 1 is non-trivial. 
In other words, we look for 5 vectors such that = Mab% admits a solution. If we now 
take Ua to be the null vector of the transpose of Mab, it is obviously orthogonal to all the 
Cq since UaCa = CaMabQb = 0. We define the non-trivial current J"^ to be: 



J' 



JI-4 + JI + JI + 24. 



(A.15) 



However, we will presently show that this current does not obey the Wess-Zumino consis- 
tency condition, therefore it cannot appear in the anomaly. 
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Wess Zumino consistency condition and J'^ 

The goal of this section is to figure out whether all possible terms that we found above sat- 
isfy the Wess-Zumino consistency conditions. To this end, we shall compute the quantities 

Oa = 5i j (fxVhu;2ViJi- §2 j d^xVhuJiViJi (A.16) 

= j d^x52{VhJi)diUJi- j d'^x6i{VhJi)diOJ2 (A.17) 

for each a = 1, .., 6. The main idea of this analysis is to find all possible linear combinations 
of the ri's such that 

6 

Y^Ca^a = (A.18) 

a=l 

If the vector space spanned by the vectors {c} is six dimensional, all J^'s are Wess-Zumino- 
consistent. If, on the other hand, this vector space is five-dimensional then we must con- 
clude that one of the J^'s is inconsistent. Since we already know that five currents can 
be generated by varying appropriate local scale invariant terms, these are manifestly con- 
sistent. Therefore the inconsistent current, if present, must be the non-trivial current of 
equation (A. 15). 

The way we shall carry out this computation is by first computing the first term in 
(A.17). The second term in (A.17) is then obtained from the first one by replacing the 
derivatives that act on oji for derivatives that act on uj2 by means of partial integration. 

We shall start with J^i. The first term in (A.17) is^^ 

S2{^ J'i) diUJi = \^{-d'uj2NR-d'AL02N) diui (A.19) 
The second term is then 

6i {Vh Ji) diUJ2 = Vh {-d'uji NR - d'Aoji N) diuj2 (A.20) 

= \/h{-d'uj2NR-V'V^{dji02N)) diioi (A.21) 

= Vh{- d'uj2NR- d'Auj2N (A.22) 

- Aw2 d'N - d' {djuj2 d^N)) diUJi (A.23) 

so that 

Qi = j d'^xVh {Au}2d'N + d'{djU2d^N)) diui (A.24) 

Similarly, from Jg: 

(52 (\//^ 4) diUJi = Vh {d'uj2 NR - Aw2 d'N) daoi (A.25) 
5i 4) diuj2 = Vh {d'uji NR - Awi d'N) 8002 

= Vh{d'uj2NR + d'{dji02d^N)) diui (A.26) 

= - j d^xVh {Au;2d'N + d'{djUJ2d^N)) diUJi (A.27) 

^^For notational clarity, notice that the variation differs by a factor of two compared to before. For 
instance, hij — >■ e"hij rather than hij — e^^hij. 
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From Jg: 



82 {Vh 4) d^oji = Vh {d'uj2 djN N + 2djUj2 d'N N) diUi (A.28) 

5i {Vh 4) a,W2 = Vh {d'LOi djN N + 2djioi d'N d^N) diL02 

= Vh {d'L02 djN dm + 2djUJ2 N d'N) daoi (A.29) 

O3 = (A.30) 



From J^: 



§2 [Vh 4) diUJi = Vh {d'uj2 AN + 2djU}2 ^d'Nd^N + Aw2 d'N) diUJi (A.31) 

<5i {Vh Ji) dii02 = Vh {d'ioi AN + 2djL0i j^d'NcPN + Aui d'N) diU2 

= Vh {d'uj2 AN + 2djL02 ^d^Nd'N - d' {djU2 d^ N)) daoi (A.32) 

^4 = y (fxVh {Auj2d'N + d'{djUJ2d^N)) diOJi (A.33) 



From J5: 



52{Vh 4) diOJi = Vh{djUJ2V'dm + d'uj2jfd^NdjN + Vjd'uj2dm) diOJi (A.34) 
61 {VhJi) diUJ2 = Vh {djcoi V'd^N + d'coi ^ff'NdjN + Vjd'ui d^ N) diUJ2 

= Vh [djUJ2 V'd^N + d'uj2 j^d^NdjN - Vj {d^'ioi d^^N)^ diUJi (A.35) 

05 = J (fxVh {Au}2d'N + d'{djU}2d^N)) diOJi (A.36) 

From Jq. 

h {Vh 4) diLOi = Vhd' {2djuj2 d^N + Awa N) daoi 

= Vh {d'AL02 N + Alo2 d'N + 2d' {djL02 dm)) dmoi (A.37) 
5i{VhJi)diUj2 = Vhd'{2dj00idm + AuiN) dii02 

= Vh {d'Au)2 N - Auj2 d'N) diUJi (A.38) 

= 2 j (fxVh {Auj2d'N + d'{djUj2d^N)) diUJi (A.39) 

We thus find that each 0^ is a multiple of 

J (fxVh {Au2 d'N + d' {djUJ2 dm) ) diOJi , (A.40) 

which means that there is one linear combination that does not satisfy the Wess-Zumino 
consistency conditions. In other words, all but one of the six J*'s can be made consistent. 
Since we have already found that five of the six J*'s can be canceled by variations of local 
terms, the one that cannot be canceled (which we called J"^) must be inconsistent. We can 
make this more precise by noticing that the consistency equation 

C1-C2 + C4 + C5 + 2c6 = (A.41) 
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describes a five-dimensional hypersurface of consistent linear combinations CqJ*. The set 
of all such Ca-vectors can be defined as those that are orthogonal to the inconsistent vector, 
Va say, such that CaVa = 0. The inconsistent vector is 

1-10112) (A.42) 



As a consistency check on our computations, notice that this is precisely the five-dimensional 
hypersurface that we mentioned above, which may be defined as all vectors that are or- 
thogonal to Ua (as defined in (A. 15)). Namely, the vector Ua is the same as the inconsistent 
vector, i.e. Ua = Va- The fact that ^7* does not satisfy the Wess-Zumino condition means 
that it cannot appear as the variation of either local or non-local terms. The fact that 
there are precisely five total-derivative terms in the anomaly, all of which can be canceled 
by variations of local terms, was also noted in [13]. 

Note: In an earlier version of this paper, we had claimed that there would be three 
instead of two possible contributions to the Lifshitz scaling anomaly, one of which was 
said to be proportional to J VhViJ^^. This erroneous conclusion has been corrected in the 
present version of this paper. 

B Role of the massive vector on the field theory side 

In this section we explore the conformal invariance of the field theory Lifshitz model from a 
different perspective. In particular, we will show that a preferred timelike vector plays 
a very similar role to the vector field appearing in the bulk. 

Our set-up is the following three-dimensional scalar model with critical exponent z = 2 

[12], 

S = [ (fxdtC = \ f <fxdt(<j)'^-{A(l)f). (B.l) 



2 „ 

The operator A is the spatial Laplacian A = 5^Widj and the dot denotes differentiation with 
respect to (imaginary) time, (j) = dtcj). The Noether current density (Ja)^ corresponding to 
the infinitesimal diffeomorphism x"" x'^ + e'^ is given via the usual definition^^ 

5eS = j d^xdtiJafdbe", (B.2) 

The current (Jf)" generates time reparametrizations and (Jj)" generates the spatial ones; 
their components are given by 

{J^y = d'<pA4>-<pd'A4> (B.4) 
(Ji)* = -<P^^<P (B.5) 
{JiY = 6{/: + did^(t>A^-di4>d^A(f> (B.6) 



*We use the notation a:* = t, i.e. the index a runs over a = t,l, 2. 



- 21 - 



where 9* = S^-'dj. One thing we see here is that (Jt)* = where £ is the Hamilto- 
nian/energy density. The conservation law reads 

dbiJa)' = -aa0(0 + AV) ~ (B.7) 

The symbol ~ denotes weak equality, i.e. equality up to terms that vanish on shell. The 
'gauge' parameter that generates the Lifshitz scaling is = 2et and e* = ex* {e is just a 
small real number). The condition for scale invariance is 

2{JtY + iJ^y = ^^{-2^'cj)A4>) (B.8) 

whose right-hand side is not zero but a total divergence. The conserved current 5" associ- 
ated to scale invariance of the theory is 

5* = 2t{jty + x\j,Y = 2t{Jty + x^{jjy (b.9) 

Note that we cannot interpret the J's as comprising an energy momentum tensor, since it 
would be far from being symmetric. 

If we couple the Lifshitz model to N and hij , we can easily write down the condition for 
conformal invariance, however the relation between the bulk (with its complete metric and 
the extra gauge field) and the field theory model is rather obscure. Clearly, the bulk metric 
does not couple to the energy momentum tensor of the field theory as defined through Jt 
and Jj, since that tensor is not even symmetric. 

So we will now make a more precise proposal about the relation between the two. We 
introduce a three-dimensional metric g^i, and a unit timelike vector so that nan"" = — 1. 
Define the projector hj' = gj' + Uafi^ , which is orthogonal to n'^, and 

= da{h^''di,(l)) + \''^h^''da(t>di,Kd (B.IO) 

then we can couple the Lifshitz model to hab and via the covariant action 

S = j (fx dt^ [{rfda(t)f - (M?) ■ (B.ll) 
This action is conformally invariant under 

5na = 2U} Ha, 5hab = 2u hab (B.12) 

This is why it is useful to introduce hJ', since the three-dimensional metric itself would 
transform as 5gab = —^uj nanh + 2uj hab (using the completeness relation gab = —nanf, + hab)- 
Of course, all of this is not very profound. We have merely replaced the spatial metric hij 
by the projection of the metric in the plane perpendicular to unit normal na- 

The claim is that (B.ll) describes the coupling of the Lifshitz model to a metric and 
a gauge field in exactly the same way as one would expect from the bulk description. 

With this fully covariant action, we can define a symmetric "stress tensor" by varying 
it with respect to gab- Due to the presence of n^, this stess-tensor is not conserved though. 
The precise equation that expresses general covariance of the theory reads 

2^^:^ = (D^Ub)^ - D, (n'^^) , (B.13) 

dgab OUb V '^^6/ 
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where on the left hand side we recognize the covariant derivative of the stress-tensor. The 
background field is the quantity that breaks the general covariance of the theory, which 
explains why this equation has a right-hand side. 

In view of (B.12), to write the conformal anomaly in covariant variables, we also need 
a variation in terms of Ua 

6S SS 

^/^A= {-4:nanh + 2hab)j + 2na — . (B.14) 

OQab OUa 

This is exactly the same as the bulk equation with Ua playing the role of the asympotic 
value of Aa- When we choose nt = N, gu = —N'^, gu = and gij = hij, the conformal 
anomaly becomes the expression we have been using all along. 

C Alternative computation of C2 

In this section we provide an alternative computation of C2 ■ Since the structure multiplying 
C2 contains R, we can take N = 1 and assume that hij does not depend on t but does 
depend on rr*. With these assumptions, the u integral separates out, yielding a factor of 
■y/vr/e. What is left is to study the operator exp(— eV^). Now we can roughly think of the 
standard heat kernel expansion as the Laplace transform of the spectral density. So by 
taking the inverse Laplace transform we can reconstruct the spectral density. The inverse 
transform of e'^ is s~^~°' /T{—a). Next, we can integrate this against exp(— es^) to obtain 



e'^/2r(-a/2) 



(C.l) 



2r(-a) 

This suggests that if the operator V has heat kernel expansion 

e"in (C.2) 



n>-l 



then has the expansion 



V- n/2 r(-n/2) 
^ ' 2r(-n) 

n>-l ^ ' 

The term with n = 1, which would contribute to the anomaly, vanishes due to the Gamma 
function. Therefore the coefficient of the B? term vanishes^^. We conclude that the coeffi- 
cient C2 in (2.29) vanishes as well. 

Notice that while this method is simpler than a direct computation, it is not powerful 
enough to determine the total derivative terms. 

D Divergent terms in the heat-kernel expansion 

The methods of section 2 allow us to compute the divergent part of the heat-kernel expan- 
sion for the model considered in equation (2.2). We present the result here for completeness: 

K(e, /, D) ~ iao(/, D) + ^~ai{f, D) + 0(6°), (D.l) 



^See also [22] for a rigorous proof of this statement. 
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where 




dt(fxNVhf{t, x){R-—AN 



(D.2) 



(D.3) 



References 

[1] D. T. Son, "Toward an ads/cold atoms correspondence: A geometric realization of the 
schrodinger symmetry", Phys.Rev. D78 (2008) 046003, 0804.3972. 

[2] K. Balasubramanian and J. McGrecvy, "Gravity duals for non-rclativistic efts" , 
Phys.Rev.Lett. 101 (2008) 061601, 0804.4053. 

[3] S. Kachru, X. Liu, and M. Mulligan, "Gravity duals of lifshitz-like fixed points", Phys. Rev. 
D 78 (2008) 106005, 0808.1725. 

[4] P. Koroteev and M. Libanov, "On existence of self-tuning solutions in static braneworlds 
without singularities", JHEP 0802 (2008) 104, 0712.1136. 

[5] M. Guica, K. Skenderis, M. Taylor, and B. C. van Rees, "Holography for schrodinger 
backgrounds", JHEP 1102 (2011) 056, 1008.1991. 

[6] S. F. Ross, "Holography for asymptotically locally lifshitz spacetimes". Class. Quant. Grav. 28 
(2011) 215019, 1107.4451, * Temporary entry *. 

[7] M. Baggio, J. de Boer, and K. Holsheimer, "Hamilton-jacobi renormalization for lifshitz 
spacetime", 1107.5562. 

[8] R. B. Mann and R. McNees, "Holographic renormalization for asymptotically lifshitz 
spacetimes", JHEP 1110 (2011) 129, 1107.5792. 

[9] M. Henningson and K. Skenderis, "The holographic weyl anomaly", JHEP 9807 (1998) 023, 
liep-th/9806087. 

[10] I. Adam, I. V. Melnikov, and S. Theisen, "A non-relativistic weyl anomaly", JHEP 0909 
(2009) 130, 0907.2156. 

[11] P. R. S. Gomes and M. Gomes, "On ward identities in lifshitz-like field theories", 1112.3887. 

[12] E. Ardonne, P. Fendley, and E. Fradkin, "Topological order and conformal quantum critical 
points", Annals Phys. 310 (2004) 493-551, cond-inat/0311466. 

[13] T. Griffin, P. Horava, and C. M. Melby- Thompson, "Gonformal lifshitz gravity from 
holography", 1112.5660. 

[14] G. Grinstein, "Anisotropic sine-gordon model and infinite-order phase transitions in three 
dimensions", Phys. Rev. B 23 May (1981) 4615-4630. 

[15] P. Di Francesco, P. Mathieu, and D. Senechal, "Conformal field theory", 1997. 

[16] D. V. Vassilevich, "Heat kernel expansion: User's manual", Phys.Rept. 388 (2003) 279-360, 
hep-th/0306138. 

[17] A. Ceresole, P. Pizzochero, and P. van Nieuwenhuizen, "The curved space trace, chiral and 
einstein anomalies from path integrals, using flat space plane waves", Phys.Rev. D39 (1989) 



1567. 



-24- 



[18] J. dc Boer, E. P. Vcrlindc, and H. L. Vcrlindc, "On the holographic renormahzation group", 
JHEP 0008 (2000) 003, hep-th/9912012. 

[19] D. Martelh and W. Mueck, "Holographic renormahzation and ward identities with the 
hamilton-jacobi method", Nud.Phys. B654 (2003) 248-276, hep-th/0205061. 

[20] S. F. Ross and O. Saremi, "Holographic stress tensor for non-relativistic theories", JHEP 
0909 (2009) 009, 0907.1846. 

[21] A. Adams, A. Maloney, A. Sinha, and S. E. Vazquez, "1/n effects in non-relativistic 
gauge-gravity duality" , JHEP 0903 (2009) 097, 0812.0166. 

[22] P. B. Gilkey, "The spectral geometry of the higher order laplacian", Duke Mathematical 
Journal 47 (1980), no. 3, 511-528. 



-25- 



